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Abstract 

A reanalysis is made for the helium abundance determination for the Izotov-Thuan 
(2004) spectroscopic sample of extragalactic H II regions. We find that the effect of under- 



lying stellar absorption of the He I lines, which is more important for metal poor systems, 
affects significantly the inferred primordial helium abundance Y p obtained in the zero metal- 
licity limit and the slope of linear extrapolation, dY/dZ. This brings Y p from 0.234 ±0.004 to 
0.250±0.004 and dY/dZ = 4.7±1.0 to 1.1±1.4. Conservatively, this indicates the importance 
of the proper understanding of underlying stellar absorption for accurate determinations of 



the primordial helium abundance to the error of SY p ~ 0.002 — 0.004. 



1. Introduction 

Izotov and Thuan (2004; hereinafter IT04) presented the primordial helium abundance 
Y p = 0.242 ± 0.002, consistently with their earlier publications (Izotov and Thuan 1998, 
and references therein), from helium recombination lines in metal poor extragalactic HII 
regions. This is significantly higher than the earlier values given by a number of authors 
(Pagel et al. 1992; Olive et al. 1997; Peimbert et al. 2000), yet is significantly lower by three 
standard deviations than the expectation from the baryon abundance constrained from CMB 
temperature anisotropics (Spergel et al. 2003) with the aid of Big-Bang nucleosynthesis 0. 
Particularly intriguing is the small errors which are common to nearly all analyses. We 
suspect that this might not represent properly the error including systematics. 

IT04 provided high quality spectroscopic data for 33 HII regions detailed enough for us 
to repeat the analysis. In this paper we are particularly concerned with the effect induced by 



: We adopt here the neutron lifetime from the conventional decay counting of ultra cold neutrons (Eidclman 
et al. 2004). 
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underlying stellar absorption on nebula helium emission lines that is poorly constrained. Our 
analysis is along the line of Peimbert et al. (2002) and Olive and Skillman (2004; hereinafter 
OS04), who introduced the helium absorption strength as a free parameter, although ours 
are somewhat more conservative. OS04 and Peimbert et al. (2000; 2002) presented an 
analysis in which plasma temperature is also determined by the helium emission lines alone. 
However, this, in principle proper approach given accurate data, induces large errors in the 
resulting helium abundance with the present accuracy of the available data for distant HII 
regions, and a trend that might be brought about with the inclusion of stellar absorption is 
buried in the noise. Since we do not see a compelling reason that plasma temperatures from 
helium and oxygen are significantly different, we take the approach that plasma temperature 
is determined by the ratio of oxygen emission lines, as was done in most of the work including 
IT04, and study the effect of stellar absorption by introducing it as a free parameter. We 
show that Y p and dY/dZ are very sensitive to the introduction of the underlying stellar 
absorption of the helium lines. 



2. Data and procedures of the analysis 

We consider 30 of the 33 H II regions given in IT04@ discarding 3 H II regions (UM133, 
Mrkl063, HS0111+2115), for which He I A4026 line is not detected. We also add NGC346 
(Region A) studied by Peimbert et al. (2000). The 31 H II regions we studied are shown 
in Table 1, where our final results are also presented. The plasma temperature of OIII is 
determined from the ratio of [O III] emissions A4363 versus AA4959,5007. The extinction 
and stellar absorption of H I Balmer lines are derived from Ha, H/3, H7 and H<5 using the 
recombination calculation of Hummer & Storey (1987) for the intrinsic line intensity ratios 
I(A)//(H/3), as 

i_m = f(\) w(R\)+a m w(m mca0 

J(H/3) F(Rf3) W(R/3) + a m W(RX) ' 1 ' 

where F[X) is the observed line intensity, W(X) is the equivalent width, ani is the stellar 
absorption and /(A)ch/3 is the extinction relative to H./3. The inclusion of H9 and higher 
Balmer lines, where available, does not modify the result beyond our interest. H8 is blended 
with He IA3889, and equation ([1]) is used to deblend the helium line. We discard H7 that is 
deblended with [Ne III]. We use the extinction curve of O'Donnel (1994), but the use of the 
different extinction curve leads only to a small difference that can be ignored here. 



2 We do not include the sample of Izotov & Thuan (1998) because the equivalent widths are not provided 
for the recombination lines. 
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We consider six He I lines, five lines for orthohelium A3890 (3s — > 2p), A4026 (5c? — > 2p), 
A4471 {Ad -> 2p), A5875 {3d 2p), A7065 (3s -> 2p), and one line for parahelium A6678 
(3d — > 2p). The A3890 and A7065 are sensitive to a fluorescent correction and thus to the 
radiative transfer. A4026 is generally weak, susceptible largely to stellar absorption. We use 
the effective recombination coefficients of Benjamin et al. (1999), which include collisional 
excitation, and the radiative transfer calculation of Benjamin et al. (2002) for fluorescent 
corrections, fx, which is controlled by the optical depth r(3890). We calculate the abundance 
of singly ionised helium as 



This is the same as that adopted by OS04 and that by IT04 up to the inclusion of eiHei- We 
do not know what ratios are to be taken for an e i for different lines. We assume here that 
all absorption strengths are identical in equivalent strengths, as was done in OS04. This 
is probably not too bad an approximation in view of the observation for absorption in B 
stars (Lennon et al. 1993; Lyubimkov et al. 2000) and the size of the resulting errors in 
our calculation that amount to £s50% of the central values: more detailed line ratios are 
the matter at a higher order level. We refer to Olive & Skillman (2001) for more detailed 
discussion for this issue. 

Our method of analysis differs from the usual one to find the parameters. We find the 
best likelihood solution in full four parameter space y + , n e (electron density), an e i and r for 
six lines, rather than minimising the sum of x 2 for each line, which has been adopted in the 
literature, for the purpose to take the correlation among the lines into account. The electron 
density is poorly constrained, but its indeterminacy affects the resulting Y little due to its 
very weak dependence in the recombination coefficient. 

We add the abundance of doubly ionised helium y ++ to calculate Y when He II A4686 
is detected. We take the oxygen abundance from IT04. The final results are displayed for 
the mass fraction of helium Y as a function of the oxygen abundance relative to hydrogen 
O/H, and linear extrapolation is employed to obtain the primordial helium abundance Y p in 
the zero oxygen abundance limit (Peimbert & Torres- Peimbert 1974), 




F(HeIA) J(Hg W(Rp) WjRelX) + a Uc i 
F{R(3) J(HeIA) W (H/3) + a m W(RelX) 



1Q/( A ) c H/3 

fx 



(2) 




(3) 



We adopt dY/d{0/E) = 18.2dY/dZ from IT04. 
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3. Results 

Extinction and stellar absorption of the hydrogen lines are determined from the 4 Balmer 
lines. The extinction parameter ch/? agrees very well with those of IT04 in spite of the 
different extinction curves used. Curiously, the stellar absorption we derived is generally 
somewhat larger than that of IT04, although we do not claim that the two are inconsistent 
since errors are generally large. For most of the case absorption strengths lie in the range of 
— 6A, but for two it is larger than 10A. The mean is am = 2.5±3.5A, where the error stands 
for the standard deviation. This is a natural value for HII regions. As noted by OS04, we 
find negative values for absorptions for 11 HII regions. This occurrence is expected from a 
large scatter in am, but might be ascribed to the collisional excitation for hydrogen emission 
lines that are not taken into account. In any case, the error induced by resetting the negative 
absorption width to zero will be small for the helium emissivity analysis, since the reference 
H/3 has a large equivalent width > 200A. 

To confirm that our results agree with those of IT04 when He I stellar absorption is 
ignored, we first carry out the analysis assuming an e i = 0. The helium abundance plotted as 
a function of the oxygen abundance [O/H] presented in Figure 1(a) confirms the trend seen in 
IT04, although individual data scatter more with our analysis. We obtain Y p = 0.234 ±0.003 
consistent with 0.2385 ± 0.0015 from the extrapolation of the results tabulated in IT04 [see 
Figure 2(a) jH Our x 2 1S rather poor: the mean is x 2 — 5.9. There are five cases (HS 
0122+0743, Mrk724, POX36, HS0128+2832 and Mrk 1236) that give x 2 > 10, implying the 
inadequacy of the procedure we assumed. 

When we include ami as a free parameter, and carry out a similar analysis, x 2 is 
enormously improved: we find x 2 = 2.0, which is dominated by a few bad fits, notably 
by one system giving x 2 — 6.1 with the zero He I absorption width. For all cases, except 
one, that showed poor x 2 with ami = 0, fits are greatly improved. In particular for the five 
cases of poor fit we noted above x 2 drops from 29.5 to 1.2 for HS0122+0743, 10.1 to 2.5 for 
Mrk724, 10.1 to 1.0 for POX36. and 18.5 to 0.9 for HS0128+2832, and from 10.6 to 4.4 for 
Mrk 1236, indicating the need to take the stellar absorption into account. 

Figure 1(b) shows the plot that can be compared with Figure 1(a) but with a non- 
vanishing ami- It is interesting to see that the derived helium abundance is more strongly 
affected for metal poor HII regions, pushing up the helium abundance extrapolated to the 
zero metallicity quite a bit and resulting in Y p = 0.250 ± 0.004. Note that the large scatter 



3 IT04 give Y for individual HII regions but do not give Y p for their new data alone. They combined with 
their 1998 data and present Y p = 0.2429 ± 0.0009. 
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around the extrapolation line seen in panel (a) is not visible in this plot any more. The \ 2 
curves are given in Figure 2(a) for Y p and for dY/d(0/H) (b) which is discussed below. 

The absorption widths derived from our fit have generally large errors, ~ 50%. The 
mean absorption equivalent widths are and = 0.40 ±0.31. There are 6 HII regions for which 
the central values of an e i are negative. The four among six are the cases for which the 
fits without stellar absorption gave unusually good \ 2 - The only cases that are marginally 
acceptable are with HS1028+3843 {\ 2 = 3.6) and Mrk35 {\ 2 = 6.1). (In the analysis with 
OHei = 0, 10 systems have x 2 greater than 6.1.) 

OS04 presented their analysis for seven HII regions common in our sample. A compari- 
son shows that our Y is always consistent with theirs within one standard deviation, although 
the errors are significantly larger in OS04 due to larger uncertainties in the temperature of 
the plasma. We find other parameters, such as temperature, are also consistent. 

Figure 3 shows the average absorption equivalent width that varies with the heavy 
element abundance. Although individual data show a rather large scatter, the trend is 
clear in this binned plot: the stellar absorption effect is more important for metal poor HII 
regions^. Since the oxygen abundance shows a tight anticorrelation with temperature of the 
plasma, this figure is also interpreted as showing the correlation of the stellar absorption with 
the plasma temperature: Higher the temperature, more the important absorption. This is 
the systematic trend that largely modifies the extrapolation of the helium abundance to the 
zero metallicity. 

In this connection another interesting quantity is dY/d(0/H), i.e., the increment of 
helium per heavy element production. The IT04 results give dY/d(0/H) = 82 ± 15 ( or 
dY/dZ = 4.5 ± 0.8, which is consistent with the final result IT04 quoted, dY/dZ = 3.7 ± 
1.2) from a sample of seven HII regions. Our analysis with an e i set equal to zero gives 
dY/d(0/H) = 86 ± 18 or dY/dZ = 4.7 ± 1.0. On the other hand, with stellar absorption we 
obtain dY/d(0/H) = 20 ± 25 or dY/dZ = 1.1 ± 1.4, which is a drastic decrease. 

We remark that IT04 use only 3 lines in deriving the helium abundance, dropping A3889, 
A4026 and A7065, with the anticipation that the other 3 bright lines are less affected by stellar 
absorption, while they use A3889 and A7065 to constrain other parameters. We also tried to 
drop the three lines when we calculate Y, but the results remain unchanged from the full 6 
line analysis. 



4 While our derived parameters are consistent with those of OS04, the effect is not discernible in OS04 
due to larger errors (introduced by the temperature uncertainty) and their smaller data set. 
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4. Conclusions 

Using the IT04 sample we showed that the neglect of the stellar absorption on the 
helium emission lines causes a large systematic effect on Y p and dY/dZ. This is due to 
the fact that the stellar absorption equivalent width shows a trend that increases towards 
metal poor systems. The inclusion of helium stellar absorption improves enormously the 
acceptance of fit to the He I emission data, especially when the quality was bad without the 
inclusion of stellar absorption. The resulting magnitude of the absorption equivalent width 
of He I A4471 line (and also of the H I Balmer lines) is on the order that is expected in a 
population synthesis calculation (Gonzalez Delgado et al. 1999) for the nebula phase. 

With the inclusion of stellar absorption we obtained the primordial helium abundance 
increased from y p = 0.234 ± 0.004 to 0.250 ± 0.004. We do not claim that the latter is 
the true value, but it is much preferred to the former on the ground of much smaller x 2 
for the fit to individual HII regions, while the 6-line fit without stellar absorption is barely 
acceptable. Or, most conservatively, one can claim that we cannot obtain the primordial 
helium abundance to the error of SY P w 0.004 or less unless underlying stellar absorption is 
properly understood. 

In our analysis we noted that the minimisation of the sum of the \ 2 over individual 
lines, ignoring the correlation among lines, is likely to underestimate the error. 

If we accept our higher helium abundance, we find ns/n-y = 7.9l 2 '4 x 10~ 10 with the aid 
of the standard Big Bang nucleosynthesis calculation (e.g., Olive et al. 2001). This baryon 
abundance is consistent with that inferred from cosmic microwave background anisotropics 
(Spergel et al. 2003). 

We also found that dY/dZ is largely affected upon the inclusion of stellar absorption: 
the original value of dY/dZ « 4 — 5 decreases to 1 ± 1. This smaller number is consistent 
with the derivative inferred form the standard solar model (Bahcall et al. 2001) AY/AZ = 
(initial — Yp) /initial = 1-4, and also with 2.1 ± 0.4 from the dwarf star atmosphere (Jimenez 
et al. 2003) using the method introduced by Pagel and Portinari (1998) who presented 3 ± 2 
for this derivative. 

This work is supported in part by Grants in Aid of the Ministry of Education of Japan 
at Kashiwa. MF received support from the Monell Foundation at Princeton. 
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Table 1. Physical Properties of HII Regions and Helium Abundance 
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HS 0735+3512 

HS 2359+1659 

HS 1311+3628 

HS 1213+3636 
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5.08 
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18.5 


+ 
-0 
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.0053 


10.6 


+ 
-0 
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.0018 


4.04 


+ 
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0062 
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7.61 


+ 
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0064 
.0039 


4.03 


+ 
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Table 1 — Continued 



Hllreion 0/H[xlCT 4 ] T e [10 4 X] C H(3 a(HI) [A] n £ [cm" 3 ] r <x(HeI) [A] Y X 2 F (w/o abs) X 2 (w/o abs) 



UM 396 


1.708 ±0.064 


1.140 ±0.015 


0.258 ± 0.028 


6.31 ± 1.22 


620 


+••■ 

620 


1.20 


+ 1.71 
-0.75 


0.00 


+0.21 


0.2471 


+0.0094 
-0.0078 


1.51 


0.2471 


+0.0081 
-0.0078 


1.52 


Mrk 1315 


1.774 ± 0.041 


1.101 ±0.008 


0.149 ± 0.017 


0.00 ±0.31 


310 


+249 
-231 


0.32 


+0.66 
-0.32 


0.15 


+0.24 
-0.15 


0.2554 


+0.0041 
-0.0043 


2.85 


0.2533 


+0.0031 
-0.0030 


3.34 


Mrk 1329 


1.787 ± 0.044 


1.079 ± 0.009 


0.178 ± 0.023 


0.82 ± 0.97 


500 


+ ■•■ 
-460 


0.48 


+ 1.24 
-0.48 


0.31 


+0.22 
-0.23 


0.2574 


+0.0055 
-0.0073 


1.71 


0.2497 


+0.0029 
-0.0023 


3.90 


Mrk 35 


1.977 ±0.058 


1.016 ±0.012 


0.230 ±0.018 


0.00 ±0.07 


335 


+270 
-237 


2.24 


+0.72 
-0.71 


0.00 


+0.08 


0.2554 


+0.0027 
-0.0023 


6.10 


0.2555 


+0.0021 
-0.0024 


6.10 



Note. - 



The dots in upward errors of n e means that one sigma value is not reached within the limit (n e = 1000) of our fit we set. 
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Fig. 1. — Helium mass fraction Y vs oxygen abundance 0/H for 30 HII regions in IT04 
(filled circle) and NGC346A (open circle) (a) without and (b) with steller absorption. The 
solid line represents the linear fit. 



-12- 




0.22 0.23 0.24 0.25 0.26 




-40 -20 20 40 60 80 100 120 



dY/d(0/H) 

Fig. 2. — Ax 2 = X 2 ~ Xmin f° r ( a ) an d (b)dY/d(0/H). The solid and dashed curves 
denote the results with and without absorption, respectively. The dotted curve shows the 
helium abundances given in IT04. A% 2 for 68% and 95% confidence levels are indicated. 
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Fig. 3. — Stellar absorption for Hel emission lines vs the oxygen abundance (filled square). 
The binned average is shown by solid circles. 



